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Abstract 

We study the neutralino decay in the supersymmetric extra U(l) models which can solve 
the /i-problem. In these models the neutralino sector is extended at least into six compo- 
nents by an extra U(l) gaugino and a superpartner of a singlet Higgs. Focussing on its two 
lower mass eigenstates x\ an d Xi, decay processes such as a tree level three body decay 
X2 ~^ Xiff an d a one-loop radiative decay X2 ~ > Xil are estimated. We investigate the 
condition under which the radiative decay becomes the dominant mode and also numer- 
ically search for such parameter regions. In this analysis we take account of the abelian 
gaugino kinetic term mixing. We suggest that the gaugino mass relation Myy ~ My may 
not be necessary for the radiative decay dominance in the extra U(l) models. 
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1 Introduction 



Recently the standard model (SM) has been confirmed in the incredible accuracy through 
the precise measurements at LEP. Nevertheless, it has still not been considered as the 
fundamental theory of particle physics and physics beyond the SM is eagerly explored. 
Along this line the supersymmetrization of the SM is now considered as the most promising 
extension However, even in this minimal supersymmetric standard model (MSSM) 
there remain some theoretically unsatisfactory features in addition to the existence of too 
many parameters. The famous one is known as the /i-problem [H. The MSSM has a 
supersymmetric Higgs mixing term \1H\H1. To cause an appropriate radiative symmetry 
breaking at the weak scale 0, we should put /i ~ 0{G F ) by hand, where Gp is a 
Fermi constant. Although in the supersymmetric models its typical scale is generally 
characterized by the supersymmetry breaking scale Ms which is usually taken as 1 TeV 
region, there is no reason why \i should be such a scale because it is usually considered to 
be irrelevant to the supersymmetry breaking. The reasonable way to answer this issue is 
to consider the origin of //-scale as some result of the supersymmetry breaking ||]. One of 
such solutions is the introduction of a singlet field S and replace /iifii/2 by a Yukawa type 
coupling XSH 1 H 2 . If S gets a vacuum expectation value(VEV) of order 1 TeV as a result 

— 1/2 

of renormalization effects on the soft supersymmetry breaking parameters, \i ~ 0{G F ) 
will be realized dynamically as /i = X(S). As is well-known, such a scenario can be 
available by introducing a nS 3 term into the superpotential and a lot of works have been 
done on this type of models f| , where the superpotential of S is composed of the terms 
XSH1H2 + kS 3 . At the price of the introduction of a new parameter k, a kS 3 term 
can prohibit the appearance of a massless axion and also guarantee the stability of the 
potential for the scalar component of S. The introduction of an extra U{l)x symmetry 
which is broken by a SM singlet field S can effectively play the same role of the introduction 
of the kS 3 term ||. A D-term for this U(l)x induces a quartic term of S in the scalar 
potential. The axion is eaten by this extra U(l)x gauge boson and disappears from the 
physical spectrum. Moreover, this extra U(l)x automatically forbids the appearance of 
LiHiH 2 in the original Lagrangian and also if we assume the unification of gauge coupling 
constants, we need no new parameter like k. Thus the models extended with an extra 
U(l)x symmetry can be considered as one of the most simple and promising extensions 
of the MSSM. Their phenomenological aspects have also been studied by various authors 
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i a is- 

The extra U(l)x models have an another interesting aspect if they are supersym- 
metrized. Their supersymmetrization introduces the extra neutralino candidates besides 
the ones of the MSSM, that is, an extra U(l)x gaugino Ax and a superpartner S of the 
singlet Higgs S. The confirmation of the extra gauge structure is one of the main parts 
of the study of extension of the SM. It is well-known that the extra U(l)x gauge struc- 
ture is often induced from more fundamental theory like superstring ||. However, recent 
precise measurements at LEP and the direct search at Tevatron suggest that the lower 
bound of the extra neutral gauge boson is rather large and it may be difficult to find its 



existence directly in near future [10|. If the supersymmetry is the true story in nature, 



there may be a new possibility to find its existence in the completely different way |l I 



Even if the mass of extra neutral gauge boson is too large to observe its existence in the 
near future collider experiment, its superpartner sector may open the window to find its 
existence. The study of neutralino sector is interesting from the view point not only of 
the investigation of supersymmetry but also of the search for the extra gauge structure. 
In particular, we should note that the gauge coupling of this extra U(l)x to ordinary 
matter fields is rather large compared with the ordinary Yukawa couplings (instead of top 
Yukawa)[] and then the neutralino sector can be substantially affected by this inclusion in 
a suitable parameter region. 

In this paper we treat the neutralino decay in the extra U(l)x models since it may be 
one of the important subjects along the above mentioned direction. The lightest neutralino 
is a candidate of the lightest supersymmetric particle. Thus if R-parity is conserved, the 
neutralino decay modes such as x° — > Xiff and X2 ~^ Xil are expected to appear as the 
subprocess of the decay of supersymmetric particles, where X2 ari d Xi are the two lower 
neutralino mass eigenstates. These decay processes have been calculated in the case of 
the MSSM under the suitable conditions |L2] . 

Recently, some attentions have been attracted to this process in relation to the CDF 
ee77 + Ifix event fll3f . Especially, related to this type of events, it seems to be a very 
interesting subject under what condition X2 ~^ Xil can become the dominant mode 
IB], III]. This is because it can give us the fruitful information on the parameters of 



1 It should also be noted that the Yukawa coupling A of XSH1H2 can be large enough compared with 
the ordinary Yukawa couplings. 
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supersymmetric models as stressed in Since this type of process is a typical one which 
may be observed in near future, its detailed study in the /z-problem solvable extra U(l)x 
models will be useful. The estimation of the widths r(x° ~ > Xil) an d r(%2 — ¥ Xiff) 
in the extra U(l)x models can be modified from that in the MSSM because there are 
new components Ax and S contained in the neutralino mass eigenstate Xi- Additionally, 
in the multi U(l)s models the abelian kinetic term mixing can occur as suggested in 
refs. IfLSl, |16| , [171 ■ As a result of this abelian kinetic term mixing, there is some changes 



in the interactions between neutralinos and ordinary matter fields [lT|. This should be 
taken into account in the analysis of these processes. Due to these effects the X2 ~~ * X?7 
dominant condition is also expected to be altered from the MSSM one. If we take the 
lesson brought from the study of the CDF type event seriously, this analysis may give us 
an important information for the model building on the additional gauge structure and 
also the Planck scale physics. 

The organization of this paper is the following. In section 2, we present the examples 
of the /i-problem solvable extra U(l)x models derived from the superstring inspired E 6 
models. After that we give a brief review of the abelian gaugino mixing whose effect 
is taken into account in the later analysis. We also examine the neutral gauge boson 
and Higgs sector to constrain the parameters of the models in terms of their present 
experimental mass bounds. In section 3, mass eigenstates and their couplings to the 
matter fields of the extended neutralino sector are studied. Based on these preparations 
the decay widths r(x° — ► Xiff) an d r(x° — * Xil) are estimated. We also study under 
what condition the radiative decay mode becomes the dominant one, which is crucially 
relevant to the CDF type event. In section 4, these decay widths are numerically estimated 
and we show what kind of parameter region is crucial for the radiative decay dominance. 
Section 5 is devoted to summary. 



2 Extra U(l)x models 
2.1 /i-problem solvable models 

There can be a lot of low energy extra U(l)x models. In these models we are especially 
interested in /i-problem solvable extra U(l)x models. From such a point of view, it 
seems to be natural to examine the models which satisfy a condition mentioned in the 
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Table 1 The charge assignment of extra U(l)s which are derived from Eq. These charges 
are normalized as Qj = 20. 

ie27 

introduction. That is, the extra U(l)x symmetry should be broken by the VEV of the 
SM singlet S which has a coupling to the ordinary Higgs doublets Hi and H2 such as 
XSH1H2. In these models the /z-scale is naturally related to the mass of the extra U(l)x 
boson and then they seem to be very interesting from the phenomenological viewpoint 
too.0 So we confine our attention to this class of models derived from the superstring 
inspired Eq models. 

There are two classes of extra U(l)x models derived from superstring inspired E 6 
models. The rank six models have two extra U(l)s besides the SM gauge structure. They 
can be expressed as the appropriate linear combinations of U(l)^ and U(l) x whose charge 
assignments for 27 of Eq are given in Table 1. There is also a rank five model called 77- 
model. Its charge assignment is also listed in Table 1. As seen from this table, there 
2 There is also a possibility that the /z-term is realized by a nonrenormalizable term X(S S / M^) 71 S H1H2 
because of some discrete symmetry In such a case (S) should be large in order to realize the 

appropriate /x-scale. As a result there is not the low energy extra gauge symmetry which can be relevant 
to the present experimental front. Because of this reason we do not consider this possibility. 
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is a SM singlet S which has the coupling XSH1H2. 77-model crearly satisfies the above 
mentioned condition. On the other hand, in the rank six models this condition imposes 
rather severe constraint on the extra U(l)x in the low energy region. In this type of models 
a right-handed sneutrino N also has to get the VEV to break the gauge symmetry into the 
SM one. If we try to explain the smallness of the neutrino mass in this context, N should 
get the sufficiently large VEV. In fact, in the case that N has a conjugate chiral partner 
N, a sector of (N, N) has a D-flat direction and then they can get a large VEV without 



breaking supersymmetry [19|. This VEV can induce the large right-handed Majorana 
neutrino mass through the nonrenormalizable term (N N) n / M^ 1-3 in the superpotential 
and then the seesaw mechanism is applicable to yield the small neutrino mass |6|, |2f| . 
However, this usually breaks the direct relation between the /i-scale and the mass of the 
extra neutral gauge boson because the VEV of N also contributes to the latter. In order 
to escape this situation and obtain the extra U(l)x satisfying our condition, we need to 



construct a U(l)x by taking a linear combination of U(l)^ and U(l) x || |16, |2(J. As such 
examples, we can construct two low energy extra U(l)x models. They are shown in Table 
2. The difference between them is the overall sign.Q In these models the right-handed 
sneutrinos have no charge of this low energy extra U(l)x- This is a different situation from 
the rank five ?7-model. Thus using the D-flat direction of another extra U(l), the right- 
handed sneutrino gets the large VEV which breaks this extra U(l) symmetry and also 
can induce the large Majorana masses for the right-handed neutrinos. This mechanism 
may also be related to the inflation of universe and the baryogenesis as discussed in ||22|| . 
As a result of this symmetry breaking at the intermediate scale, only one extra U(l)x 
remains as the low energy symmetry. We will concentrate ourselves on these three U(l)x 
models (X = 77, £±) in the following study. 

We focus our attention to the minimally extended part of these models with an extra 
U(l)x and a SM singlet Higgs S. Other extra matter fields like color triplet fields (g,g) 
and the right-handed neutrino N, which are introduced associated with the extension, 
are irrelevant to the present purpose and we can neglect them. Thus the relevant parts 
3 As discussed in refs. || ^o), can also be obtained only by changing the field assignments for Q x . 
This insight allows us to construct new models, which can induce an interesting neutrino mass matrix 
PH by using the charge assignments Q x and for the different generations |2(J]. However, in this 
paper we shall not consider such models for simplicity. 
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Table 2 The charge assignment of the extra U(l) x which remains unbroken after the VEV 
of N becomes nonzero. They are obtained as Q^ ± = dh ^ ± tQx- 



of the superpotential and soft supersymmetry breaking terms are 

W = XSH 1 H 2 + h u QU c H 2 + h D QD c H 1 + h E LE c H 1 + ---, 
C soft = -J2m^ + (AXSH 1 H 2 + h.c+---) 

i 

1 3 . . 

+ -(M w J2*w*w + M Y*YXY + M x \ x \ x + M YX \ Y \ x + h.c), (1) 
1 0=1 

where <pi represents the scalar component of each chiral superfield contained in the models. 
M W ,M Y and M x are the gaugino masses.f] We assume the Yukawa coupling A and soft 
supersymmetry breaking parameters to be real, for simplicity. 

2.2 Abelian gauginos mixing 

Next we briefly review a particular feature in the neutralino sector caused by the abelian 
gauge kinetic term mixing in the supersymmetric multi U(l)s models. In supersymmetric 
4 In this expression we introduced the abelian gaugino mass mixing as Myx , which might exist as the 
tree level term at the Planck scale and also be yielded through the quantum effects. 
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models gauge fields are extended to vector superfields 

v wz (x, e, e) = -ea^ev** + mex - mex + \qqqqd, (2) 

where we used the Wess-Zumino gauge. A gauge field strength is included in the chiral 
superfield constructed from Vwz in the well-known procedure, 



W Q (x,0) = (DD)D Q V WZ 
Here we should note that W a of the abelian gauge group is gauge invariant itself. In terms 



UX a - A6 a D + m^a,aUd^V v - d»V») - A99a ua0 d^. (3) 



of these superfields the supersymmetric gauge invariant Lagrangian can be written as 

£ = 1 (W a W a ) F + (V ex P (2g QV wz )^) n , (4) 



32 

where $ = ((f), if), F) is the chiral superfield and represents matter fields. Its generalization 
to the multi U(l)s case is straightforward. The supersymmetric gauge kinetic parts are 
obtained by using chiral superfields W° and for U(l) a x Z7 (l)b as 

— (w aa W*) + — (W ha w£) + *^(W aa W*) . (5) 
32 V >f 32 V J f 16 V a /F v 1 

Here we introduced the mixing term between the different Z7(l)s. This can be canonically 

diagonarized by using the transformation, 

/ 





(6) 



1 — tanx 
^0 1/ cosx 

This transformation affects not only the gauge field sector but also the sector of gauginos 
A a 5 and auxiliary fields -D a ,&-Q As easily seen from the form of the last term in eq. 
(4), the change induced in the interactions of gauginos with other fields through this 
transformation can be summarized as 

9aQaX a + g° b Q b X b = g a Q a X a + (gabQa + 9bQb) X\ (7) 

where X a ^ are canonically normalized gauginos. The charges of U(l) a and U(l)b are 
represented by Q a and Qb- The couplings g a , g a b and gb are related to the original ones g® 
and g® as, 

9a = 9a, 9ab = ~9° tan X, 9b ~ " 



cosx 



5 This shift in the D-term changes the scalar potential and can affect the symmetry breaking at the 
weak scale. However, we will not refer to this problem here. 
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These coupling constants at weak scale will be determined by using the renormalization 

. However, such 



group equations from the initial values at the high energy scale jig , | 
a study is beyond our present purpose and we will treat them as the parameters in the 
later analysis. 



2.3 Neutral gauge sector 

In the previously introduced extra U{l)x models, the gauge symmetry of the electroweak 
sector at the low energy region is SU{2)i J x U(1)y x U(1)x- In order to bring the correct 
symmetry breaking for these models, we assume that Higgs fields get the VEVs as follows, 



(Hi) 



where v\ + v\ = (246 GeV) 5 







[') 
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\ 1 




\v 2 j 



(S) = u, (9) 

\ V 2 J 

; 2 ) is assumed. For simplicity, all VEVs are assumed 



to be real. Under these settings in order to constrain the parameters of the models, we 
investigate some features of the gauge boson sector. 

For this purpose we need to determine the physical states at and below the weak 



scale ||17|| . The mass mixing between two neutral gauge fields appears associated with the 
spontaneous symmetry breaking due to the VEVs of eq. (9) around the weak scale. In 
the present models the charged gauge sector is the same as that of the MSSM. In the 
neutral gauge sector we introduce the Weinberg angle 9w in the usual way,Q 

= cos 6 w Wl - sin 9 W B^ = sin 8 W W^ + cos 9 W B^. (10) 

Here we used the canonically normalized basis (Z^, X^) so that has already been 
decoupled from [Z^X^). The mass matrix of the neutral gauge fields (Z^X^) can be 
written as 
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Am 2 
cosx'' 



m x = m z s w tan \ + 2 Am sw 
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M 2 Z , 
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'121 



In the following we use the abbreviated notation sw = sin 0w and cw = cos 



W ■ 
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In this expression m%, Am 2 and M§, represent the values of corresponding components 
in case of no kinetic term mixing (x = 0), 

m z = \(9w + gl)v 2 , 

Am 2 = l -{g 2 w + g 2 Y ) 1/2 9xv 2 (Qi cos 2 (3 - Q 2 sin 2 (3) , 

Ml = \g\ (Q\v\ + Q\v\ + Q 2 s u 2 ) . (13) 

The mass matrix eq. (11) can be diagonalized by introducing a mixing angle £. The 
canonically normalized mass eigenstates are written by using \ and £. Their concrete 
expressions are given in Appendix A. We also present there the interaction Lagrangian of 
the neutral gauge bosons and matter fermions for the later use. 

The mixing angle introduced for the diagonalization of the mass matrix eq. (11) is 
given by 

-2cosx(m|svysinx + Am 2 ) 

tan 24 = —5 — — - : 5—5 — 5 - — . (14) 

M-2,1 + 2Am z sw sin x + m z s w sm X ~ m z cos X 

In general the mixing angle £ is severely constrained to be small enough by the precise 
measurements at LEP [jTU] . From the study of radiative symmetry breaking it has been 



known that tan/? ^ 1 is generally favored. In fact, it has been shown in ref. that the 
suitable radiative symmetry breaking could occur for 1.4 ^ tan/3 < 2.1 in £_ model. We 
will adopt 

tan/? ~ 1.5 (15) 

as its typical value throughout this paper. Therefore, in case of sin \ = 0, since Qi cos 2 (3 ~ 
Q2 sin 2 (3 is not satisfied in the present three models, we need to consider the possibility 
that the small £ is realized because of Am 2 <C M|/, which is equivalent to f 2 ,f| u 2 . 
If sin x 7^ 0, however, there may be a new possibility to satisfy the smallness of £ even if 
Am 2 Mf / is not satisfied. Such a situation can be expected to occur if the following 
condition 

sinx 2~~ = 1 2 i 9X 2\M2 — (<5i cos 2 /? - Q 2 sin 2 /?) (16) 

is valid. In this case Q 1 cos 2 /3 ~ Q2sin 2 /5 is not required unlike the sinx = case but 
instead of that the tuning of sinx becomes necessary. The constraint on the value of u 
also becomes very weak. Since this possibility for the small vn^ x compared with m| is 
interesting enough for the explanation of the smallness of |£|, we will also consider the case 
with such a mixing angle sin x in the following discussion as one of the typical examples. 
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The present model-independent bound on the mixing angle £ is |£| < 0.01 p4|. If 
we impose this bound on the models, we can restrict the allowed u range in each model. 
Here it should be noted that the mixing angle £ has no A dependence. In order to show 
this constraint coming from the neutral gauge sector, we plot the contours of the mixing 
angle |£| = 0.01 for each model in the (sinx, \u\) plane in Fig. 1. The lower regions of 
the contours are forbidden in each model. It is noticeable that rather small value of \u\ is 
generally allowed in ^-models in comparison with 77-model. From this figure we find that 
the kinetic term mixing sin% can affect the lower bound of \u\ substantially. In 77-model 
the larger sinx reduces the required bound of \u\ values. In ^-models 7 there are special 
values of sin% which make the lower bound of \u\ very small, as anticipated in eq. (16). 
Thus in these models the rather light extra Z° may be possible.0 

Related to the fact that rather large \u\ is generally required except for the case with 
the special sin x value, it will be useful to remind again the origin of /i-scale in the present 
models. In these models the vacuum expectation value u is relevant to the /z-scale. Based 
on this feature we may need to put the upper bound on A to keep \i to be the suitable 
scale from the viewpoint of radiative symmetry breaking as discussed in |J. If we use 
the present Higgs mass bounds, however, A can be effectively constrained as shown in the 
following subsection. 

2.4 Higgs sector 

Higgs sector is changed from that of the MSSM due to the existence of the singlet S 
and its coupling XSH 1 H 2 to the Higgs doublets Hi and H 2 . Its brief study can give us 
some useful informations of the allowed region of the parameter space|25, ^6[. If we take 



account of the abelian gauge kinetic term mixing, the scalar potential for the Higgs sector 
can be written as, 

, 2 



V = \(<& + &){& + 

7 £±-models have the symmetric feature mutually with respect to the sign of sinx so that they are 
expected to show the similar behavior in their phenomenology. This comes from their characteristics of 
the charge assignments. 

8 In this case the extra U(l)x gaugino is also expected to affect largely the rare phenomena like 
(1 — > e7 and the EDM of an electron fill]. 
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if , / 2 2\ , 9X 

+ - <g Y tanx{v 1 - v 2 ) + 
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cos x 

+X 2 vlvl + X 2 u 2 vl + X 2 u 2 v 2 + m^ 2 + m\v\ + m|w 2 - 2AAwit> 2 , (17) 

where Qi, Q2 and Qs represent the extra U{l)x charges of Higgs chiral superfields H\, H 2 
and S. At the minimum of this potential, the mass matrices for the Higgs sector are given 
as follows respectively, 
• charged Higgs scalar sector: 
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neutral Higgs scalar sector: 

/ m|cos 2 /3(l + 4) + AA«tan/3 ^ sin2/3(-l + il^+^l) _ AXu 
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(19) 



where g = Jg^ + g Y and we define 0, C2 and (3 as 



> , . 9xQi , 
Ci = ^ytanxH , C2 



, . 9xQ2 
-9y tan xH , 



C 3 



(20) 



cos x cos x cos x 

The overall factor of a mass matrix of the charged Higgs sector is somehow changed 
from that of the MSSM due to the coupling XSH1H2. However, the mass eigenstate of 
charged Higgs scalar can be obtained in the same form as the MSSM case 



= sin (3H± + cos (3H^* 



(21) 



and its mass eigenvalues are expressed as 



M 2 H± = m 2 z c 2 w 1 



2A 2 



+ - 



AXu 



gly J sin (3 cos (3 



(22) 



The A 2 term is added to the MSSM one and then the charged Higgs mass takes smaller 
value than that of the MSSM for the same value of fi = Xu. On the other hand, the 
neutral Higgs mass matrix is too complex to be diagonalized analytically. However, if 



9 It may be useful to note that the sign of Ci> C2 and £3 is reversed between £ + -model with sinx and 
£_-model with — sinx. 
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we note that the smallest eigenvalue of the matrix is smaller than the smallest diagonal 
component, we can find the tree level upper bound of the lightest neutral Higgs mass. By 
diagonalizing the 2x2 submatrix at the upper left corner of eq. (19) we can obtain] 



2 ^ 2 

m h0 < m z 



2A 1 

cos 2 2(3 + -= T sin 2 2(3 + -= T (& cos 2 (3 + C2 sin 2 (3)' 

9w + 9y 9w + 9y K 



(23) 



The first two terms correspond to the bound which is derived from the usually studied 
model extended with a gauge singlet. 

As easily seen from these results, these Higgs masses have the crucial dependence on 
A and u. One of the important difference between the present models and the MSSM 
comes from the fact that the //-term is replaced by the Yukawa coupling \SH\H2- If we 
impose the present experimental bounds on the Higgs masses, useful constraints can be 
obtained in the (A, u) plane. The present mass bounds on both of the charged Higgs and 
the lightest neutral Higgs are ~44 GeV [|27|. We use this bound and show the allowed 
region in the (A, u) plane in Fig. 2. Since it is found to be insensitive to the models 
and also the sinx value, we take £_-model with sinx = as an example. Here for the 
lightest neutral Higgs we used the result obtained by the numerical diagonalization of the 
mass matrix eq. (19). It should be noted that only the u > region is allowed. This is 
completely dependent on our choice (A > 0) for the sign of 

Additional important constraints on fi can be obtained from the condition in the 
(fi,Mw) plane coming from the search of the neutralinos and charginos at LEP ||28|| . If 
we assume tan/3 ~ 1.5, the allowed region in this plane is roughly estimated as,[^] 

Id, M w > 40 GeV (for Am > 0), 

11 (24) 
\fjt\, M w > 100 GeV (for \u < 0). 

The chargino sector in the present model is not altered from the MSSM and then these 

conditions on \i can be used as the constraint for A and u. Thus the allowed region of the 

(A, u) plane are found to be determined by the lower bound of the lightest neutral Higgs 

10 The A and u dependence of the Higgs mass eigenvalues is included in the terms, which are composed 

of Au and even powers of each of them. Thus the sign of u is related to that of A. Here it should also be 

noted that in the present notation u > corresponds to the ordinary /i < case. 

n It should be noted that this restriction has been derived under some assumptions, for example, the 

gaugino unification relation My = § tan 2 OwMw- However, we will apply them for the general My and 

Mw here. These constraints correspond to the condition for the chargino mass m\ ^ 65GeV. 
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mass for all models. It corresponds to the surrounded region by the dashed lines in Fig. 
2. If we combine this with the result obtained from Fig. 1, we can restrict the allowed 
region in the (A, u) plane for each model with a certain sinx value. We will use this fact 
later. 



3 The decay width of X2 i n ^° Xi 
3.1 Neutralino sector 

In this subsection we examine the structure of the neutralino sector and also define the 
mass eigenstates of the charginos and squarks/sleptons sector, which are necessary for 
the calculation of the neutralino decay. Starting from the superpotential and soft su- 
persymmetry breaking terms given in eq. (1) and using the canonically normalized 
basis defined by eq. (6), we can write down the modified quantities from the MSSM, 
which are relevant to the neutralino sector, that is, the neutralino mass matrix and the 
gaugino-fermion-sfermion interaction terms. If we take the canonically normalized gaug- 
ino basis M T = (—iXy V ,—iXY,—iXx,Hi,H2,S) and define the neutralino mass term as 
£™ ass = -iAfTMAf+ h.c, the 6 x 6 neutralino mass matrix M. can be expressed as 

( M w m z cw cos (3 —m z cw sin (3 ^ 



(25) 



My C\ —mzsw cos (3 mzsw sin (3 

C\ C2 C3 C4 C5 

mzcw cos (3 —mzsw cos (3 C3 Am Afsin/3 

—mzcw sin (3 mzsw sin (3 C4 \u Xv cos (3 

^0 C 5 Xv sin P Xvcos(3 J 

where v and u are defined by eq. (9). Matrix elements C\ ~ C5 are components which 
are affected by the kinetic term mixing. They are represented as 

n w f Myx n a# f 2 , M x 2M YX sinx 

d = -M Y tanx H , C 2 = M Y tan x H 5 5 , 

cos x cos x cos X 

C 3 = 4= \9y tanx + ) vcosp, C 4 = ( ~9y tanx + \ us in/5, 

V2 V cosx/ V2 V cosx/ 

n 1 9xQs , na , 

Cs = 75^7" (26) 
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Neutralino mass eigenstates xl (i — 1 ~ 6) are related to A/} through the mixing matrix 



U<ij dbS 



X° = U T J\f . 



(27) 



The change in the gaugino interactions can be confined into the extra U(l)x gaugino 
sector and new interaction terms can be expressed as, 



ip* f -g Y Y tan x + 9x ® x \ 
V cos x 



i 

+H* (-g Y Y tan x + 9x 9 x ) \ X H - ( -g Y Y tan X + XxHH 



-g Y Y tan x + 9x ® X j \ x iptp 
cosx J 



cosx 



cosx 



(28) 



where ifj and ip represent quarks/leptons and squarks/sleptons. Higgs fields (H\, H 2 , S) 
are summarized as H and the corresponding Higgsinos (Hi,H 2 ,S) are denoted as H. 
The charges of U{1) Y and U(l)x are denoted as Y and Q x . As a result, the parts 
corresponding to the gaugino component of the neutralino x°-fermion-sfermion vertices 
are represented by the following factors 



SliYQx) 
Z?(Y,Qx) = - 



1 

V2 
1 

V2 



gwU U T 3 + g Y YU 2i + ( -g Y Y tanx + 9x ® x ) U 3i 



cosx 



g Y YU 2i + ( -g Y Y tan X + ] [fo 



(29) 



where the suffixes L and i? stand for the chirality of the coupled matter fields ip and their 
charges are defined in terms of the left-handed chiral basis as presented in Tables 1 and 
2. 

Additionally, it is also useful to define the chargino and squark mass eigenstates here 
for the forthcoming calculation. Taking account of eq. (1), the chargino mass terms are 
given as 

/ 

^mass = ~~ (-^1 ) — ^ 



—\u \/2mzCw cos (3 

\f2m z c w sin (3 M w 



(30) 



The mass eigenstates xt are defined in terms of the weak interaction eigenstates through 
the unitary transformations, 



(xt 

•0+ 

\ X2 



I 



w(+)t 



-iA H 



v X2 



( 



Hi 



(31) 
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(32) 



S quarks and sleptons are also relevant to the neutralino decay. When we consider 
this subject, all flavors can be treated in the same way except for the top sector. If they 
appear in the internal lines, the stop may be especially important because of the largeness 
of its Yukawa couplings and then we only consider the stop sector in such cases. However, 
in the neutralino decay modes which contain the ordinary fermions in the final states, top 
quark is too heavy to be included in them and it is irrelevant to such processes. 

In the following analysis we do not consider the flavor mixing in the squark and slepton 
sector, for simplicity. Thus the sfermion mass matrices can be reduced into the 2x2 form 
for each flavor. This 2x2 sfermion mass matrix can be written in terms of the basis 
(II, Ir) as 

\m f \ 2 + Ml + Dl m f (A f + XuRf) 
m* f (A* f + XuR f ) \m f \ 2 + Ml + Dl 

where rrif and M\ R are the masses of ordinary fermion / and its superpartners /l,r, 
respectively. We assume M\ R is universal for all flavors. Rf is cot (3 for up-sector and 
tan/? for down-sector. Soft supersymmetry breaking parameters Af are the dimensionful 
coefficients of three scalar partners of the corresponding Yukawa couplings. D\ and D\ 
represent the D-term contributions, which are modified in the present models as follows 

D 2 L = ± X -m\ cos 2/3(1 - (1 ± Y)s 2 w ) + \g\Q'AQ>X + Q^l + QW), 

D 2 R = -\m 2 z s 2 w Y cos 2(3 + -^iQ'xiQ'A + Q' 2 v\ + Q' s u 2 ), (33) 

where the upper sign in D L corresponds to the up-sector sfermions and the lower one to 
down-sector sfermions. The primed charge Q' x stands for the modified charge due to the 
kinetic term mixing and defined as gxQ'x = ~ 9yY taxix + 9xQx/ cosx- We should note 
that these D-term contributions cannot be neglected in the extra U(l) x where u tends 
to be large. In such cases it will be useful to note that the positivity condition of the 
sfermion masses may induce no condition on the soft scalar masses. We define the mass 
eigenstates (/i,/^) as 





(34) 

Under our assumption for the reality of soft SUSY parameters the above chargino and 
sfermion mass matrices are real and then and V? become the orthogonal matrices. 

By now we have finished the preparations for the calculation of neutralino decay in 
the present models. If R-parity is conserved and the lightest neutralino is the lightest 
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superparticle, the decay of the next-to-lightest neutralino X2 i n t° the lightest neutralino 
x\ can be expected to appear in the various superparticle decay processes. As the rep- 
resentative decay modes of x\ Xi> the tree level three body decay \\ — > Xiff an d 
the one-loop radiative decay X2 ~^ Xil have been calculated in the MSSM framework 



12| , |T4] , |29|| . In these studies, which decay mode of these becomes dominant has been 
shown to be crucially dependent on the composition of X2 an d Xi an d then on the SUSY 
parameters. It is very interesting that one-loop decay mode can easily dominate the tree 
level process in the suitable parameter region. As was recently stressed in ref. [O], if 
the CDF type event relevant to x[> ~~ > X?7 happen to be observed dominantly instead of 
X2 ~^ Xiff> it can constrain the SUSY parameter space severely. In the following part of 
this section we shall analyze the decay widths of x\ ~ * Xiff an d X2 ~^ Xil i n the present 
extra U(l)x models and also qualitatively discuss the condition on the SUSY parameters 
for the X2 ~^ Xil dominance. 

There exist other decay modes like two body decay into the lightest Higgs X2 ~^ 
and the cascade decay mediated through the chargino as X2 ~^ Xi( e ^e) — ► XiC^e(ez/ e ). 
If the h° is light enough for the threshold to be opened satisfying m^o — m^o > m h o, 
the first one can be a relevant mode. The second one may not be suppressed if xt is 
lighter than x\ even in the case that x\ is composed of the same ingredients as the case 
where X2 ~ * Xiff is suppressed. Although these points should be taken into account 
in the analysis, through the numerical calculation of the mass eigenvalues h° at least 
seems to be heavy enough not to open the threshold in the parameter region (Xu, mz) 
where we are interested in. For the chargino mediated cascade decay the threshold can be 



opened but the existence of its suppression mechanism has been pointed out in ref. |14| . 
Therefore, in this paper we concentrate our attention on the comparison of X2 ~* Xi7 an d 
X2 ~^ Xiff- F° r this purpose we shall firstly calculate the decay width of both modes. 
We are particularly interested in the case of rather small neutralino masses since in such 
a case these neutralino decays may be observed in the experiment in near future. 



3.2 xl^xlff 

There are two types of diagrams which contribute to the tree level three body decay. 
They are shown in Fig.3. Top quark cannot be a final state so that the contribution 
from diagram (b) is generally suppressed by the small Yukawa coupling. The phase space 
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integral can be analytically done in the limit that the mass of the final state fermion / is 
zero. This seems to be generally rather good approximation and we adopt this result of 
the phase space integral in the present estimation. Thus the decay width for this process 
can be expressed as^ 



96vr 3 



4 4 
m j ~ m i ( 4 , 4 

— - — ^ — mZ + m\ 



8m?m^ + 24m,-mf In — - 
l > J rrii 



(35) 



a=l 



where the vertex factors T a can be expressed by using the mixing matrix element Uij in 
the neutralino sector as 
1 



Am 2 Zi 



9 {1) + 



9x 



4m| 2 



9 {2) + 



2 cosx 



9x 
2 cosx 



Qisin£ U±jU, 



-g 



(i> 



9x 



2 cosx 



Q 2 sin£ ] U 5j U 5i 



9x 



2 cosx 



Qs sm£U 6j U 6i 



Qi cos£ U Aj U Ai 



-9 {2) + 



9x 
2 cosx 



, F (i) 



-Qs cos ^U ej U ej 



/ii 



(/=tf) 



E 



2 cosx 

g x )c/ 5 , + zj(y; Qx)c/ 5 i) 

+ Q x )Zi(Y, Q x ) + z£(y, Qx)Zg(y, Q x )) V/ a V 2 f 



(36) 



where 



,(2) 



9w_ 

2cw 
9w 

2c 



w 



(cos £ + svk tan x sin £) , 
(— sin£ + svk tan x cos £) . 



(37) 



and JF 2 comes from the diagram (a). The mass eigenvalues of the neutral gauge bosons 
are expressed as and mz 2 - Mj is mass eigenvalues of sfermion mass matrix eq. (32). 

Jo 

The effective neutral current couplings etc. are deviated from ones of the MSSM due 
to the existence of the extra U(l)x and the abelian gauge kinetic term mixing. Their 
concrete expressions are presented in Appendix A. Z^Y, Qx) and Z^(Y, Qx) are defined 
by eq. (29). The diagram (b) gives T{± and !F^ D ' E ' is obtained by replacing Z^(Y, Qx), 
~Z*(Y,Q X ) and U 5j in T { J =U) with Z^(Y,Q X ), Z^(Y,Q X ) and U 4j , respectively. 

12 It should be noted that in the limit of mr — > there is no interference term like Ft Ft between 

J JL JR 

the different fermion chiralities in T\ (a = 1,2). 
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It is useful to examine under what condition this decay width can be suppressed based 
on eqs. (35) and (36). As was noticed up to now ||14|| , there can happen the dynamical 
suppression depending on the composition of and Xi which is determined by the SUSY 
parameters. For the contribution from T\ and T<i they are suppressed unless both of X2 
and Xi are dominated by the Higgsinos. In JF 3 and T± there are the contributions from 
both of the Higgsino and gaugino components in x\ an d Xi an d then it seems to be difficult 
to expect the suppression due to the neutralino composition. However, there is the crucial 
suppression due to the small Yukawa coupling and also the small left-right mixing V( 2 in 
the sfermion mass matrix. These features can be summarized as follows. There appears 
the dynamical suppression effectively in such a case that one of x\ an d Xi is dominated 
by gauginos and the other is dominated by Higgsinos. Although this is the same as the 
MSSM situation, there is a noticeable feature in the present extra U(l)x models. In case 
of the S dominated neutralino, it has no mixings with Ayp and Ay. Moreover, it has no 
couplings with ordinary fermions. If this is the case, it is not necessary for the gaugino 
dominated neutralino to be an almost pure photino in order to suppress this three body 
decay unlike the MSSM. Later this point will be discussed in more detail again. 



3.3 x°2^xh 

Next we proceed to the calculation of the one-loop radiative decay X2 ~~ * Xil- This 
has already been studied in the MSSM framework [I2|. From the gauge invariance, as 
suggested in P0 |, it is easily found that the effective interaction describing this process is 
given as 

£eff = G^XiF^. (38) 
Using this effective coupling Q, the decay width is written as 

where and rrij are the masses of Xi an d X°j- Our main problem is the estimation 
of the effective coupling Q. One- loop diagrams contributing this coupling are given in 
Fig. 4. In the diagrams (la) and (lb), only the stop contribution cannot be neglected 
because of its large Yukawa coupling. After some algebraic manipulation, it is obvious 
that this coupling can be obtained as the coefficient of $ ■ / terms where and e M are 
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the momentum and the polarization vector of photon. In case of small neutralino masses 
m it j <C m,M where m and M respectively represent masses of fermions and bosons in the 
internal lines, the neutralino mass dependence disappears from these one-loop amplitudes. 
Its only dependence on the neutralino sector comes through the mixing matrix of the 
neutralino sector. The effective coupling Q can be summarized as follows, 



32tt 2 



M: 



E -/Hf)^+ E 



a=l,2 



4/( 



m t 



m 2 t 



a=l,2 



m c 



mi 



mt 



w - 



+ 3J( 



G1 + ^ J( 



(40) 



///,-, ' ' ///,-, ' j ' 2m a nit 

where m a and M H ± stand for the masses of the charginos and the charged Higgs. The 
charged Higgs mass expression is presented in eq. (22). The first and second summa- 
tions should be taken for the stop mass eigenstates and the chargino mass eigenstates, 
respectively. Each term with a vertex factor Gf comes from Feynman diagram numbered 
with % in Fig. 4 and their concrete expressions are presented in Appendix B. Kinematical 
functions /(r),/(r) and J(r) are defined as 



m = 

I(r) 
J(r) 



r 
1 



1 + 



lnr 



2(1 -r) 2 
1 



2(1 



1 -r 
1 +r + 

-3 +r 



2r 
1 — r 
2 



1 



lnr 
— lnr 



(41) 
(42) 

(43) 



For checking this formula, we assume that Xu, M w , M Y , M x <C m z and the stop mass 
matrix is diagonal (V a p = 5 a p). In such a case, for the situation is the same as the 

MSSM and they can be taken as, 13 
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(-) 



w}p = = = = 0. 



Mass eigenvalues of charginos are approximately written as 

m\ = V2mzCw cos (3, m\ = V2mzCw sin (3. 



(44) 



(45) 



For Uij, if we put gx = and A — > but keeping Xu) constant, can be approxi- 
mated as 



Uu = sw, Uii = cw, Uij = sin/5, U$j = cos/9, other = 0, 



13 Here the sign conventions are taken so as to make both mass eigenvalues positive. 
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Zf{Y) = -V2g 2 s w Q em , Zf{Y) = V2g 2 s w Q em , Zf(Y) = Zf{Y) = 0. (46) 

Using these expressions, it can be easily checked that Q is reduced to the MSSM result 
calculated in this parameter setting [T2| . 

The feature of eq. (40) is rather similar to the one of the MSSM. As easily seen from 
the structure of Gf in Appendix B, there is not the special neutralino configuration in 
which the drastic suppression mechanism works for T(x° — > Xil) unlike T(x° — > Xiff)- 
This is an important feature to consider the neutralino decay processes. 



3.4 The radiative decay dominant condition 



As was clarified through the study of the CDF event ecyy + fix |13| , the neutralino 
decay can give the valuable information on the SUSY parameters. Based on the naive 
perturbative sense, as x 2 Xil is the higher order process compared with x 2 ~^ Xiffi 
the former is expected to be largely suppressed by the small couplings compared with 
the latter. However, in the present case the neutralinos are complicatedly composed 
of the various ingredients and two decay modes imply the different feature depending on 
their compositions which are determined by the SUSY parameters. If the signature of the 
radiative decay mode is dominantly observed, the SUSY parameter space can be strictly 
restricted due to the suppression condition of the tree level three body decay. Thus it will 
be useful to study how this situation can be changed in the extra U{l)x models. 

For this investigation it is convenient to rewrite the neutralino mass matrix eq. (25) 
in terms of the usual photino and Higgsino basis which is often used in the MSSM case. 
It can be written as 
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M w s 2 w + M Y c 2 w 


(Mw — My)swcw 


cwCi 
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(M w — M Y )sw c w 


M w c w + M Y s w 
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C3 cos f) — C4 sin fi 


C3 sin f) + C4 cos fi 
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C3 cos f} — C4, sin fi 


—\u sin 2/3 


\u cos 2/3 















C3 sin fi + C4 cos fi 


\u cos 2/3 


\u sin 2/3 


Xv 




V 








C 5 





Xv 





/ 



(47) 

14 It has been suggested that there is also a kinematical suppression of the three body decay when x° 
and x? are nearly degenerate m| — m| rrij fl4|| . However, in our study we will not refer to such a 
parameter region. 
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where we define the neutralino basis of this matrix as (— zAi, — 1X2, — , H a , Hb, H c ). 
Throughout this study we assume that the gaugino masses Mw and My take the smaller 
value than 200GeV. 

In the MSSM case the radiative decay dominant condition is expressed as [O, [14[ 



My, tan/? ~ 1. (48) 



The second one is natural from the viewpoint of radiative symmetry breaking and we 
assume that it is satisfied in our study as mentioned before. The first one is nontrivial 
but it may not be necessarily required strictly in some parameter region as pointed out in 



ref. As easily seen from the part of eq. (47) corresponding to the MSSM neutralino 
sector, we find that in the MSSM with the condition eq. (48) the almost pure photino Ai 
and the one of Higgsinos become the lower two neutralino mass eigenstates as far as 
M w ,M Y ,Xu <C m z . This situation realizes the suppression of the three body decay as 
discussed in the last part of subsection 3.2. On the other hand, this kind of suppression of 
three body decay seems not to be realized in the present extra U(l)x models even if the 
above condition is satisfied. This is because of the existence of the extra U(l)x gaugino 
which has the mixings with every neutralino components. Thus in order to suppress the 
tree level three body decay it is necessary to resolve this mixing effectively and produce 
the purely Higgsino-type neutralino. Although various possibilities may be considered, 
we are particularly interested in the case with M w 9^ My. 

The first possibility is to make Ai and/or A2 decouple from one of the Higgsinos by 
imposing 

Ci ~ 0, u > v, (49) 

in addition to eq. (48). The first one requires My sinx = Myx and it is always satisfied 
in case of no kinetic term mixing. The second one should be usually satisfied in the extra 
U(l)x models to overcome the small mixing condition on £ as discussed in the previous 
part. As shown in Tables 1 and 2, C3 sin /3 ± C4 cos (3 ~ cannot be satisfied in the 
present extra U(l)x models. However, if u is large enough, C5 becomes large and as a 
result C 3 sin /3±C 4 cos (3 ~ can be effectively satisfied. Under this situation the Higginos 
Hb can decouple from Ai and A2. The value of A is related to which neutralinos become 
the lower two neutralino mass eigenstates and then it seems not to be severely restricted 
by requiring the radiative decay dominace. As easily seen from the above mass matrix, 
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A3 ans H c tend to decouple from other fields under the condition (49) and the situation is 
reduced to the MSSM one. The feature of x 2 an d Xi is expected to be similar to the one 
of the MSSM. When the composition of these states is interchanged, the same suppression 
is also expected to occur. In this possiblity it should be noted that M w ~ My will not be 
necessarily required like the MSSM as far as Mx takes the similar value as Mw and My. 
In £±-model with the suitable sin% value, the large u is not necessarily needed. In such 
a case, although the Z' becomes rather light, the radiative decay dominance can not be 
expected. In this case sin \ — seems to be preferable for the radiative decay dominance. 

The second possibility is to make the lightest neutralino be the almost pure S. As 
mentioned in the sebsection 3.2, S has no mixings with Aw and Ay and also no couplings 
with ordinary fermions. Thus if we consider the situation that the next-to-lightest neu- 
tralino is the mixture of Aw and Ay and the lightest neutralino is dominated by S, the 
three body decay can be suppressed. This gives a new window which does not require 
the condition Mw ~ My. A very light neutralino dominated by S is considered in the 



different context in ref. [31]]. To realize this situation it is necessary to impose 

C 2 > C B , Ci ~ 0, Xu > m z . (50) 

The first one means that Mx needs to be rather large compared with u. We need a partic- 
ular supersymmetry breaking mechanism which can realize the large hierarchy among soft 
gaugino masses such as Mx S> My. If u ^> v which is generally the preferable situation 
for the extra U(l)x models, the next-to-lightest neutralino is almost the mixture of Ai 
and X 2 (i.e. Aw and Ay) and also the lightest neutralino H c which is purely S. Starting 
from this case, we can get other composition for the lightest neutralino which realize the 
radiative decay dominace by shifting the values of Mx and u. If we assume u > v, the 
lightest neutralino becomes the mixture of Hb and H c . This situation can be realized in 
the £-t-models with the suitable sinx value as found from Fig.l. If the condition C 2 3> C5 
is changed into C5 ^> C 2 v which is equivalent tow > Mx v, the lightest neutralino 
becomes and the situation becomes similar to the MSSM case except that x 2 does not 
need to be a photino like state but is enough to be any states composed of Aw and Ay. 
It should be noted that these new possibilities are related to the large //(> mz) and/or 
sinx ^ case, where Ax and S can play the crucial role. In the sin% 7^ case, C\ ~ 
requires the existence of nonzero Myx- The validity of this condition should be checked 
by using the RGE study in each model. 
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4 Numerical analysis 



The arguments in the previous section are qualitative one on the suppression mechanism 
for the three body decay of x° compared with the radiative decay. It is necessary to 
proceed the numerical calculations to treat the subtlety of the parameter dependences 
and also restrict in more quantitative way the SUSY parameter space where the radiative 
neutralino decay becomes the dominant mode. As suggested above, there may be a 
new window of the SUSY parameters in the present extra U(l)x models and it may be 
possible to escape the constraint eq. (48) on the gaugino mass in the MSSM. To clarify 
this we compare the two decay modes numerically. In the study of this direction the most 
interesting parameters are the gaugino masses. In addition to them, u and A will be also 
important in the present models because it is relevant to the extra Z° mass and also the 
//-scale. 

Before going to the numerical analysis of these decay widths, it will be useful to 
summarize the allowed parameter region. We have already presented the contraints on 
A and u in Figs. 1 and 2. By combining these results, for the typical values of sin% the 
allowed region of u is roughly estimated as, 
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(51) 

where Myy > 40 GeV should be satisfied. Here we should note that sinx affects the 
neutralino decay widths eqs. (35) and (39) not only directly through the vertex factors 
and the mixing matrix but also indirectly through determining the lower bound of u. For 
the soft supersymmetry breaking parameters we assume the typical values as follows, 

A — Af — 200 GeV, M L = M R = 200 GeV. (52) 
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Additionally, gy = gx and Myx = are also assumed.^ The gaugino mass Mx is 
treated as a free parameter and also the gaugino masses My/ and My are assumed to take 
not so large values such as 40GeV < M\y,My < 200GeV. Under this parameter setting, 
the branching ratio Br = T(x2 Xil) /\F(X2 ~~ * Xiff) + ^(xl ~^ Xil)] are studied in 
the (A, u) and (M W ,M Y ) planes for the typical values of sinx and M X - Through this 
study we found that the decay width of the radiative decay is in the rather wide range 
O(10~ 6 ~ 10~ 10 )GeV depending on the parameters. Although from the viewpoint of the 
experimental detectability it may be possible to restrict furthermore the parameter region 
based on the absolute value of r(x° ~~ ► Xil)-> we are interested mainly in the radiative 
decay dominance conditions and then we focus our attention only on the Br value here. It 
should also be noted that Br gives the same value for £_-model with sin% and £ + -model 
with — sin \- 

At first we examine Br under the condition of My/ = My = Mx{ ^ 200GeV) in the 
(A, it) plane. As an example, we take £_-model which has the rather small lower bound 
of u. In this model it is expected that there is no severe restriction on the value of A. 
In fact the numerical studies show that Br > 0.98 is realized almost through all the 
region which satisfies the constraints coming from Figs.l and 2, although for a certain 
\i value around ~ 600 GeV there is a shallow valley where Br gives the slightly smaller 
value compared with other region. That valley moves in the (A, u) plane by the order 
of Xu ~ O(10 lr " 2 )GeV following the change of the value of sinx from —0.2 to 0.2. This 
shift originated from the change of sin x becomes larger as Mx becoming larger. When 
Mx becomes larger, Br < 0.90 occurs at the small A region such as A ^ 0.2. These 
qualitative features are found to be common to all models. The difference between 77- 
model and £-t-model is that the latter can have the smaller bound of u. As a result, for 
the same value of A, /i in £±-models can take the smaller values than that in 77- model. In 
such a small /x region Br has the tendency to become smaller as far as the small gaugino 
masses are assumed. This is because the gaugino- Higgsino mixing cannot be extracted in 
the lower lying neutraino eigenstates. Anyway we can safely conclude that the radiative 
decay dominance is good enough in the whole region of (A, u) as far as Mw = My = Mx 
is satisfied. 

15 Although these should be determined in terms of RGEs analysis, we make these assumptions only 
for simplicity. 
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Next we proceed to the study of Mw and My dependence of Br. For this purpose 
we estimate Br in the (Mw, My) plane. In Fig. 5 we show the results for £_-model as an 
example. The global feature of this kind of plot seems to be characterized by the value 
of «(= Am) if M x is fixed. In the case of Am < m z (Figs. 5a, 5b and 5c), M w ~ M Y 
seems not to be severely required. This point has been already pointed out in the MSSM 
case [l4|. However, in this model the larger violation of the relation Mw — My seems 



to be allowed compared with the MSSM case. When Mx becomes larger compared with 
Mw and My, the Br > 0.9 region shrinks into the smaller Mw, My region and also there 
appears the new Br > 0.9 region in the large Mw, My domain, where Mw ~ My is not 
required. These behavior of Br may be understood as follows. Accompanied with the 
change of Mx, the level crossing occurs between Xi and x° and then their ingredients are 
interchanged. And in the region of Mx where the separation between Xi and X2 * s large 
enough, Br > 0.9 is realized. In the case of Am > mz (Figs. 5d, 5e and 5f), when Mx 
is smaller compared with mz, the Br > 0.9 region appears as the beltlike zone around 
the Mw ~ My line but the width of this region is not so narrow. This means that 
the next-to-lightest neutralino should be the almost photino Ai to realize the radiative 
decay dominance and then Mw — My is preferable. Under this condition the mixture 
of A 2 , A 3 , H a , H b and H c can decouple from Ai. As M x becomes larger, the Br > 0.9 
region has the tendency to occupy the wider space where Mw — My is not required. 
The reason of this Br behavior can be understood from the qualitative arguments in the 
previous section. Although we show here the results only for one model, we have checked 
that other models also showed the similar qualitative features. So these results can be 
considered as qualitatively general ones. 

Finally we would like to stress that in the extra U(l)x models there is the wide 
parameter region where the radiative decay becomes the dominant mode of the neutralino 
decay. This region contains the new possibility such that the relation Mw ~ My is 
completely violated in comparison with the corresponding parameter space to the case of 
the MSSM |14 |. This can be possible because of the existence of Xx and S. The neutralino 
decay may give us various informations on the extra gauge structure. 
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5 Summary 



We studied the decay of the next-to-lightest neutralino into the lightest neutralino in the 
extended models with an extra U(l)x and a SM singlet Higgs S, which can solve the 
/^-problem as the result of its radiative symmetry breaking. In this study we took account 
of the abelian gaugino kinetic term mixing. At first we investigated the neutral gauge 
sector and Higgs sector in order to constrain the parameter space of the models. Through 
this analysis we showed that the VEV (S) and the Yukawa coupling A of the singlet Higgs 
S were constrained in the suitable region. Next the width of the one-loop radiative decay 
and the tree level three body decay were calculated. Based on those results the suppression 
condition of the three body decay were qualitatively discussed and we suggested that there 
could be a new possibility to escape the constraint on the gaugino masses M w ~ M Y for 
the realization of such a suppression in the MSSM. This is due to the existence of the 
extra U(l)x gaugino and the singlet field S. For more quantitative analysis the branching 
ratio of the radiative decay was numerically estimated in the (A, u) and (M w , M Y ) planes. 
As a result we found that the /z-problem solvable extension with the extra U(l)x could 
largely modify the parameter space which realize the radiative decay dominance from 
that of the MSSM. Especially, it was pointed out that the condition Myy ~ My for the 
gaugino masses is not necessarily required for the radiative decay dominance as far as Mx 
is large enough. In the extra U(l)x models the slepton and squark decays which contain 
the above processes as the subprocesses can be largely affected by the existence of the 
extra gauge bosons and the singlet Higgs. These results seem to be interesting for the 
future accelerator experiments. In the supersymmetric models the extension with extra 
Z7 (l)s may have the interesting and fruitful phenomena in their superpartner sector and 
its extra gauge structure may be seen through the study of the superpartner sector. The 
further study of this aspect will be worthy enough. 
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Appendix A 

In this appendix we give the concrete expressions of the interaction Lagrangian of the 
neutral gauge sector. Original states which are not canonically normalized are represented 
by the mass eigenstates (A 1 *, Zf, Z$) as 

A» = A li - c w tan x (sin ^ + cos£Z£) , 

Z^ = (cos£ + sw tan x sin ^) + (— sin £ + sw tan x cos£) Z% , 

x» = ^izr + s^zn, (53) 

cos x cos x 

where A 1 stands for the real photon field and Zf is understood as Z^ observed at LEP. 
Using these mass eigenstates, the interaction terms of these gauge fields with ordinary 
quarks and leptons in this model can be expressed as, 

tf } = F h! JlIvJl + F%> f R ^f R , 

= F^hlM + Fjgfwfa (54) 



where the coefficients Fj l j etc. are defined as 



F fL = P ~ 2 QomS 2 w ) g {1) - eQ cmClv tan x sine + -^—Q% sin£, 
JL \ / 2cosx 

FP = -2Q cni s 2 w g {l) - eQ cm c w tan x sin f + - 9X Q f x sin £, 
1 icosx 

Ff} = (r 3 - 2Q cm s^) - eQ cm c w tan X cos f + 7 ^—Q , j cos£, 
^ v ' 2cosx 

F? } = -2<5 em s^5f (2) - eQemCw tan x cos £ + - ^ X cos £■ (55) 
J K 2 cos x 



Qx and stand for the U(l)x charges of f L and f 



R- 
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Appendix B 



We give here the concrete expressions of the vertex factors Gf (i = 1 ~ 5) in eq. (40). 



huVlv la (u 5j Z^) - U 5l Z^(^ - huVlv 2a (u bi Zf{~) - U 5j Z t R (~ 

QwW^wlfJ (U 4i U 5j - U 5l U 4j ) + (z£(-l)z£(l) - ^(-1)^(1)) 

-gwW^wtJ (y 5j z&(-i) - U 5i Z^(-l)) - gwwtJw^ (u 4i z^(i) - U 4j Z^(l) 

+Asin 2 /3 IgwWgwg (U 4j U 6i - U 6j U 4i ) + (z&(-l)U 6j - Z^{-\)U { 



, (56) 



+Acos 2 /3 [gwW^wiJ (U 5i U ej - U 6t U 5j ) + W^W^ (z^{l)U &i - Z^(l)U 6j )] , (57) 



w£^wiJ (U u U 4j - U 4i U l3 ) + w\^wf a n {U u U l3 - u u u 5j ) 

+ 7! w ™ w ™ V* u *i - u ^j) 



(+)tw(-)t 



(58) 



Gi 



9W 



cos/3 (g w W^w£ (U u U 4j - U 4j Uij) 

+^w[JwtJ (u 4j z^{-\) - U 4i Z^(-l)) - ^W^W™ (U 5j U 6i - U 6j U 5i ) 

(UijZ^i-1) - U lt Z^(-l)) - \w[t ] W^ {U^Um - U 6j U u )) 
+ sin/3 (-gwW^wg? (U lt U 5j - UmUij) 

+ ^Wia )t Wia ) (u 5j Z^(l) ~ UxZ^l)) - A VF (-)t VF (-) {UijU& . _ Uf .. U4i) 



(59) 

where in these equations we abbreviate the U(l)x charges in the expression of Z^(Y, Q x ) 
and Z- l (Y : Q x )- G% can be obtained by making the replacement such as sin/3 — > cos/3 
and cos (3 — > — sin (3 in G% ■ 
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Figure Captions 
Fig. 1 

The allowed region in the (sinx, \u\) plane due to the constraint on the mixing angle 
£ between the extra U(l)x and the ordinary Z°. The contours of |£| = 0.01 for three 
models are drawn. £ + and rj- models correspond to solid, dashed and dotted-dashed 
lines, respectively. The lower region of each contour is forbidden. 

Fig. 2 

The allowed region in the (X,u) plane for £_-model with sinx = 0. The contours of the 
present mass bounds of the lightest neutral Higgs scalar and the charged Higgs scalar are 
shown by the dashed and dashed-dotted lines, respectively. The surrounded region by 
the dashed lines and the upper region of the dashed-dotted one are allowed. The solid 
lines represent the boundary (A : \u = 40GeV and B : \u = — lOOGeV) coming from 
the experimental searches of charginos and neutralinos. The region sandwiched between 
them is forbidden. 

Fig. 3 

Diagrams contributing to the tree level three body decay x° - ► Xiff- 

Fig. 4 

One-loop diagrams contributing to x° — > x^7- The chirality flip occurs at the fermion 
internal lines and/or Yukawa vertices. In Figs. (2a) and (2b) we show the representative 
ones. 

Fig. 5a 

The contours of the branching ratio Br = 0.9,0.7 and 0.5 of £_-model with sinx = 
in the (Mw,My) plane, which are represented by solid, dashed and dashed-dotted lines, 
respectively. Parameters are set as A = 0.15, u = 600GeV and M x = 50GeV. 

Fig. 5b 

The same contours of Br as Fig.5a. Parameters are set as A = 0.15, u = 600GeV and 
M x = 400GeV. 
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Fig. 5c 

The same contours of Br as Fig.5a. Parameters are set as A = 0.15, u = 600GeV and 
M x = lOOOGeV. 

Fig. 5d 

The samw contours of Br as Fig.a. Parameters are set as A = 0.5, u = 600GeV and 
M x = 50GeV. 

Fig. 5e 

The same contours of Br as Fig. 6a. Parameters are set as A = 0.5, u = 600GeV and 
M x = 400GeV. 

Fig. 5f 

The same contours of Br as Fig. 6a. Parameters are set as A = 0.5, u = 600GeV and 
M x = lOOOGeV. 
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